Active Motion Assisted by Correlated Stochastic Torques 
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The stochastic dynamics of an active particle undergoing a constant speed and additionally driven 
by an overall fluctuating torque is investigated. The random torque forces are expressed by a 
stochastic differential equation for the angular dynamics of the particle determining the orientation 
of motion. In addition to a constant torque, the particle is supplemented by random torques which 
are modeled as an Ornstein-Uhlenbeck process with given correlation time Tc- These nonvanishing 
correlations cause a persistence of the particles' trajectories and a change of the effective spatial 
diffusion coefHcient. We discuss the mean square displacement as a function of the correlation time 
and the noise intensity and detect a nonmonotonic dependence of the effective diffusion coefficient 
with respect to both correlation time and noise strength. A maximal diffusion behavior is obtained 
if the correlated angular noise straightens the curved trajectories, interrupted by small pirouettes, 
whereby the correlated noise amplifies a straightening of the curved trajectories caused by the 
constant torque. 

PACS numbers: 05.40.-a,87.16.Uv,87.18.Tt 



I. INTRODUCTION 

Recent experimental studies evidenced that living be- 
ings, such as bacteria [ll, different water fleas fish 
[1| , birds [1] and insects [3, @| ) are able to sustain a con- 
stant mean speed over large time scales. Thus, these 
animals do not only perform passive Brownian motion, 
but instead exhibit a so-called active stochastic motion 
Q, occurring far from equilibrium. 

Such systems form a large field of recent theoretical in- 
terest [1,0, including the influence of fluctuations. One 
might relate the sources of fluctuations on animal motion 
either in the food supply [lo| or in the propulsive engine 
[ll'] that drives the unit or by assuming a molecular ag- 
itation imposed externally outside [1]. In this approach 
one widely uses similarities of the animal motion and 
Brownian motion which Einstein and Smoluchowski in- 
dependently described in a probabilistic theory [l^, [13] ■ 
Later on, Langevin encoded this theory of stochastic sys- 
tems into the Newtonian notation of equations of motion 
[l3 | , which present a first formulation of a stochastic dif- 
ferential equation. This concept, completed by an ac- 
tive element describing the self-propulsion, is nowadays 
widely used for the description of fluctuating animal mo- 
tion f7Hll|. Of foremost interest are the specific effects 
of the random impacts on the motion. 

Our focus here is on agents subjected to a torque, re- 
sulting in circular motions. This feature can be caused, 
for example, by an external magnetic field with profound 
physiological implications [lB| , or due to boundary condi- 
tions [3|. The agents themselves may exhibit a preferred 
turning direction, either caused by asymmetries in the 
propulsion as occurring in the chemotaxis of sperm cells 
[IgI \v?\ , as a search strategy [H, [l^ , or also by interac- 
tion with other agents, which leads to the emergence of 
swarming characteristics [1, [2^ . 

We start out by studying the planar motion of ac- 



tive particles with a constant velocity vq subjected to 
a constant torque $7. To take into account thermal or, 
more generally, the stochastic influence of the surround- 
ing, we first consider an additional white Gaussian noise 
^(t) driving the angle dynamics and discuss the resulting 
effective diffusion coefficient D^s (2ll - [2^ . Next, we pro- 
pose a more realistic model considering some persistence 
in the motion. Therefore, we replace the white Gaus- 
sian noise by a time-correlated Gaussian noise (colored 
noise), namely an Ornstein-Uhlenbeck process (OUP) 
9{t) |25l43C|. which implies a memory for the curvature 
feature. This generalization then allows applications 
reaching from stochastic polymer dynamics |31| over spi- 
ral wave motion 32 1 up to animal trajectories [^. This 
present study using a constant torque Q. can principally 
also be applied to a different problem, namely the con- 
trol of diffusion of electron beams in drift chambers in 
presence of additional magnetic fields We want to 

mention that other noise statistics, especially white Pois- 
sonian noise [llj, [s^ i34], are conceivable. Their study 
would form an interesting future extension of the present 
work. 

In a driven motion with a constant speed vq the cor- 
relations with correlation time Tc generate a persistence 
length Ic of the trajectory: = vqTc- As will be shown 
below, the inclusion of correlations will hence character- 
istically modify the dynamics. The resulting dynamics 
will be quantified by the spatial diffusion coefficient I?cff ■ 



II. AGENTS DRIVEN BY WHITE GAUSSIAN 
NOISE 



Our starting point is the equation of motion in two 
dimensions with a constant velocity vq and a random 
torque with the mean value fi. The dynamics for the po- 
sition vector f{t) = {x{t),y{t)} derives from its velocity 
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vector, reading 
d 



dt 



f{t) = voe^(^t) = Vq (^cos (j){t), sin0(t)^ , (1) 



with (j){t) denoting the orientation of the velocity vec- 
tor. This orientation if governed by a constant 
torque and supplemented by random fluctuations, yields 
a stochastic dynamics for (j){t) 



at vq 



(2) 



We consider the noise £,{t) to be a white Gaussian noise 
with the corresponding noise intensity denoted by and 
with vanishing mean. The 1/vq in front of the random 
force is caused by the fact that the tangential acceleration 
scales as vqc/). It expresses the circumstance that fast 
agents cannot change their direction as quickly as slower 
ones. 

For a vanishing noise intensity the constant torque 
leads to circular shaped motion, whose size is dictated 
by their cyclotron radius R — vo/n. A non- vanishing 
noise intensity on the other hand induces a more erratic 
behavior. The mean square displacement (MSD) corre- 
sponding to this dynamics can be directly calculated by 
using the shifted Gaussian transition probability distri- 
bution P(A(/), r). Then the general expression for the 
MSD of a particle with a constant velocity uo, i.e. 



it- 



(cos(A0(T)))dT, (3) 



yields in the case of white Gaussian noise 



/ white 



cos [VLt)At. 



(4) 

For small times t a ballistic behavior cx results while for 
large times a crossover to diffusive motion cx t takes place. 
The effective diffusion coefficient in two dimensions 
is related to the MSD via the well-known relation: 



lini 



At 



(5) 



As a consequence we integrate Eq. in the limit of large 
times. Furthermore, we substitute a dimensionless inte- 
gration variable x = D^Vq^t and obtain together with 
the new parameters Dq — Vq/{2D^) and luq = ^v^/D^ 
an expression for the effective diffusion coefficient, read- 



ing 



7^ white 



Do 



' COS {ujqx) dx 



Do 



Df 



f72 



(6) 
(7) 



Thus, we recover the result for a white Gaussian angle 
drive with constant torque. This finding was obtained 



before in Ref. [2i| and has been re-addressed with Ref. 
[121, see for more detailed discussions [23j]. Obviously, 
the diffusion as a function of the angular noise exhibits 
a maximum, with Z?cff = I'd/ (4^) D^ = v^fl, and 
vanishes for both D^^ — and — )■ oo. This is con- 
trasted to the behavior of the diffusion without a torque, 
Dos^ = ^^o/(2-D?) 0,[36|, which diverges for 0. 

The explanation is as follows: for small noise intensi- 
ties, agents subjected to a constant torque move in circles 
around a fairly quasi-stationary center, whereas in ab- 
sence of a torque the propagation proceeds along almost 
straight lines. With noise included, both motion types 
perform diffusion. This leads to a reduction of the mean 
square displacement in case of zero torque. In contrast, 
agents subjected to a constant torque start to spread over 
the space which is expressed by the numerator in Eq. ([T]) , 
which grows linearly in D^. For large noise intensities D^, 
this leads to a suppression of diffusion as accounted for 
in (H). 




FIG. 1: (Color online) Effective diffusion coefBcient D^^'^" 
versus intensity of the Gaussian white noise for different 
torque strengths il, with the constant velocity vq = 0.5 and 
= 0.01 within theory, see in eq. ([7|. 

We note that a finite torque always decreases D^g^^^, 
because the cosine-contribution in Eq. (g]) suppresses the 
value of the integral. This difference becomes most strik- 
ing for small noise intensity D^. In Fig. [T]the behavior 
of eq. ([7]) is depicted as a function of for several 
ri- values. 



III. AGENTS DRIVEN BY AN 
ORNSTEIN-UHLENBECK-PROCESS 

A. General outlines 

As a physically more realistic extension of our model, 
we next study time-correlated noise (i.e. colored noise) 
for the angular drive instead of delta-correlated noise 
around its mean Q. Accordingly, we use an OUP as 
the most natural continuous- valued colored noise l26l430ll . 
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Hence, our new system with the constant velocity vq and 
the torque n is described by the modified angle dynam- 
ics, reading 



e(t) - - -e{t) 

Tn 




(8) 
(9) 



This colored noise dynamics assumes one additional aux- 
iliary process, namely the OUP 6{t) with the correlation 
time Tc, the Gaussian white noise ^(t) and the corre- 
sponding noise intensity D^/t^ . 

Note that we included here a dependence of the cor- 
relation time to the power n in the noise intensity im- 
plying a corresponding change in dimension for I?j as n 
changes. In doing so, several different physical features 
can be addressed with no need to perform additional cal- 
culations. For example, for n = the random variable 
0{t) merges with the Wiener process W{t) in the limit 
of infinite correlation times Tc — >■ oo. In specific detail: 

= ^ m = ./2D-^fat')dt' = y^^wit). 

On the other hand, the case with n = 1 delivers a vari- 
ance of 6, which is independent of the correlation time 

^From a theoretical point of view, the case with n ^ 2 
is of utmost interest. The limit Tc — )■ yields Gaussian 
white noise, i.e. 9{t) — > ^y2D^£^{t), and thus allows a 
comparison with the findings taken from Sec. HH 

The corresponding ^-correlation function reads in the 
stationary limit 



-l-n 



(10) 



depicting a strong n-dependence of the variance, i.e. the 
factor in front of the exponential. 

To study the angular dynamics <j){i), we use the results 
of Ornstein and Uhlenbeck [25| , which were later general- 
ized by Chandrasekhar to motion in a higher dimensional 
space 37]. Therewith we can write the probability dis- 
tribution in the stationary limit for a transition during r 
from (pi to 4>2 and conditioned by 6i as 



X exp 



vl [(t>2 - 01 - rir - OiTc {l - e 



(11) 



Here, the mean square increment of the angle holds 



D(t) = 



2t 



3 -|- 4e '■c — e 



(12) 



This resulting Gaussian expression can be directly used 
to calculate the mean square displacement of our angle 



dynamics from Eq. ([3]) , yielding with statioary Gaussian 
distributed initial values 0i and 9\ 

((f(t) - ro)2)ouP = 2vl f\t - r)e-*M cos (r!r)dr . 
Jo 

(13) 

The function in the exponent reads 



*(t) = *c 



wherein we define 



(14) 



(15) 



as a parameter which scales monotonously with Tc and 
-D^ and can therefore be used to discuss the asymptotic 
behavior, later on. 

Without constant torque (fl — 0) and with n = this 
reproduces the result of [3^, where the motion of fish 
was analyzed within a similar model. Furthermore, Eq. 
([TS]) exhibits a crossover from a ballistic behavior oc to 
diffusive motion cx t. This can be seen by studying the 
following two temporal limits 



d 

t-^oo dt 



lim ^{{r{t) - fo)^)oup = 



(16) 

POO 

2vl / e"*(^) cos (f7T)dT. 

(17) 



The right-hand side of Eq. pT| is obviously a non- 
vanishing constant which can be identified with the effec- 
tive diffusion coefficient via Eq. ([5]), yielding our central 
finding, namely 



^off -y 



-*(t) 



COS (57T)dT. 



(18) 



For further analytical discussions it will be helpful to 
use a dimensionless representation of Eq. (|18p . To this 
end, we substitute x = ^'^^ and introduce the rescaled 
variables 



yielding 



r)OUP 



2T'i-''De 



and 



„*o(l- 



(19) 



COS {ijOcx)dx. (20) 



In comparison with the results for a white Gaussian angle 
drive [c.f. Eq. ([T])] we notice different definitions of the 
constants Dc and lJc and an additional exponential in the 
integrand, namely e*''^^^'^ ' The latter converges 
to unity for ~^ and increases monotonously with 
^'q. For n = 2, the parameters Wc and coincide with 
those defined in the white noise case, respectively, with 
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Do and loq- In this case ^'o describes the deviation from 
the results under white noise. 

Also note that the integral in Eq. (f20|) can be ex- 
pressed in a serial expansion. Substituting first z — 
^'o exp(— x/^'o) as a new variable and expanding after- 
wards the remaining exponential under the integral in a 
Taylor series, one can calculate the integral in each sum- 
mand, yielding 



oo 



nOUP _ Tj *, 



1 



kl 



fe+i 



(21) 

The asymptotic behavior of the diffusion coefficient can 
be readily derived: For vanishingly small ^E'o, the first 
item in the sum with k = dominates. Hence we obtain 



lim 

*0-!-0 1 



(22) 



which in particular recovers the white noise result of eq. 
([7]) for the case n = 2. 

The opposite asymptotics of large ^Pq can be inspected 
by looking at the original integral definition in Eq. ^U\i. 
(see also 38] ) for the case without an applied torque and 
n = 0. In this limit the major contributions to the in- 
tegral stem from values at the lower boundary x ^ 0. 
Upon expanding the double exponent until the second 
order finds 



lim 



oc lim 

'i'o—^oo 



exp - 



(23) 



which tends to zero much faster compared to the case 
without torque with ujc — 0. 



B. Asymptotic behavior in the absence of torque: 

n = 




FIG. 2: (Color online) Mean square displacement {{r{t)—f'o)'^) 
for the dynamics given with Eqs. ([8]) and ((9]) for = 0, = 
0.01, vo = 0.5 and with n = 0. Simulations (points) agree 
well with the theory (Eq. (|13|) : dashed lines). The left bar is 
oc t'^, the right bar is oc t. For Tc > 3, the ratio of subsequent 



crossover times behaves according to tcr 
^/T0, otherwise according to tcross(Tc)/fc 



cross 

(lOr, 
ss(10r,) = 10^. 




First we study the crossover between the ballistic and 
the diffusive behavior. Towards this aim we expand the 
two exponentials in Eq. (|13p. A case differentiation, 
due to the two corresponding exponential time depen- 
dencies, yields two parameter dependent regimes where 
the crossover between ballistic motion and diffusion is 
realized 

„2 n-2 

if ^-0 < 1 ^ tcross « , (24) 



/„2 «-l 

if *0 > 1 ^ tcross « W . (25) 

These crossovers are corroborated with our simulation 
results and the numerical evaluations of Eq. (|13l) in Fig. 
[21 We see that the ratio of crossover times where \l/o > 1 
holds (with the assumed parameters this is true for > 
3), behaves according to icross(Tc)/tcross(10rc) = \/T0. 
In case ^E'o < 1, the ratio of crossover times increases 
severely to icross(T'c)Across(10Tc) = 10^. 



FIG. 3: (Color online) The effective diffusion coefficients Dcff 
are depicted as a function of the correlation time Tc for several 
n values within simulations (points) and theory (Eq. (|13[) : 
lines). Other parameters as in Fig. (2] While decreasing for 
n = and approaching a non-zero constant for n = 1, DcH 
diverges cx for n — 2, see Eq. (|26| . 

The different crossover behaviors thus reflect the dif- 
ferent scaling properties of DcS for the different regimes 
of V&o in agreement with the result that 

ff *o < 1 ^ « -\f^, (26) 

,,3 

ff *o >i ^ « iv^- ^^^^ 

This all is in perfect accordance to the numerical results, 
as depicted in Fig. |31 We notice the strong qualitative 
changes in the motion of agents which are induced by 
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different n. While for n = and n = 1, the effective 
displacement decreases with the correlation time, it even 
increases with Tc for n = 2. That the diffusion coefficient 
increases unbounded with cx can be seen also in Eq. 
pop where, as mentioned in the last paragraph above, 
the integrand increases monotonously with ^E'o, which is 
the only parameter dependent on Tc for n ~ 2. 

On the other hand, regarding the limit Tc — > 0, D^s di- 
verges for n < 2. Only for n = 2, the diffusion coefficient 
remains finite and converges to the result of Mikhailov 
and Meinkohn [36,] . We notice that the effective noise in- 
tensity to which our particles are subjected (i.e., 
increases with n for correlation times smaller than one, 
and decreases with n for correlation times larger than 
one. For Tc = 1 the diffusion coefficient of the three cases 
merge to the same value. In case of n = 1, 2 for smaller 
correlation times an effectively larger noise changes the 
orientation of the particles repeatedly compared to n = 0. 
That is why the diffusion coefficient is reduced and is 
largest for n = 0. On the other hand, the effective noise 
intensity acts in the opposite way for Tc > 1 leading to a 
reduced diffusion in case of n = 0. 

Figure U] depicts for n = 1 different trajectories for a 
vanishing torque and n = I within the Tc region where 
Dctf changes only marginally (see Fig. Thus, we rec- 
ognize qualitative changes of the motion, which are just 
induced by the growth in correlation time, while the vari- 
ance of the OUP [Eq. (|10l) ] stays constant. One sees that 
the structure of the trajectory changes from an uncorre- 
lated random sequence to a waltzing like motion where 
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FIG. 4: Spatial trajectory within an angular OUP dynamics 
with constant variance n — 1 for different correlation times 
Tc while the parameters hold as in Fig. [7] but with SI — 
0, n = 1 and the time length ti = 1000. The trajectories 
changes qualitatively from a random way being unoriented 
each moment to an agent performing a waltzing-like dance. 
The effective diffusion coefficient remains unchanged for the 
shown cases (cf. Fig. |3}. 
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FIG. 5: Spatial trajectories of an OUP-driven angle without 
additional torque {Q = 0) at constant n = 2, = 0.01, and 
vq — 0.5 for different correlation times Tc and the time length 
ti — 1000. We see a slight growth in the effective displacement 
for the larger correlation times as also presented in Fig. [S] 

longer straight pieces of the trajectory are interrupted by 
left- and right-turning pirouettes. We can interpret this 
behavior as an increment of the trajectories' persistence 
as Tc increases. 

FigureOdepicts trajectories for n = 2, while = still 
holds. Conforming to the behavior shown in Fig. [3l the 
displacement slightly increases for the larger correlation 
time, due to the decreasing variance [see Eq. (ITUl) ]. which 
induces a straighter motion. Therefore, the curves also 
appear smoother for increasing Tc. 

C. Phenomena of applying finite torque: Q 

We next turn to the effects arising from the application 
of a constant torque. The effective diffusion coefficients 
are displayed in Fig. [5] for different n-values, together 
with simulation results as a function of increasing color 
Tc- We observe that the time correlation induces upon 
varying Tc a maximal mean square displacement at a fi- 
nite value of the correlation time for all three cases. A 
similar graph could be also displayed for the dependence 
on the noise intensity (not shown, see also below). 

In case that n = 0, where the integrated noise inten- 
sity in the OUP dynamics becomes independent of the 
color Tc, the maximum is most pronounced. Correspond- 
ing trajectories with constant torque are shown in Fig. [T] 
The curvature of these trajectories reflects the combined 
influence of the constant torque and of the temporal cor- 
related random torque forces inherent in 0{t). 

For small color the motion is dominantly circular 
with small random perturbations. This can be under- 
stood in view of Eq. 0, since small correlation times 
imply strong relaxation so that random changes become 
essentially white noise with 

e{t) ^ ^fw^m • (28) 

Consequently, in the cases that n — Q and n = 1 the 
white noise sources tend to vanish. The trajectories fol- 
low the motion with the constant torque f2 along trajec- 
tories with a fixed curvature cx and therefore 



FIG. 6: (Color online) Effective diffusion coefficient Dcff ver- 
sus correlation time Tc for several n- values, a constant torque 
Q. = 0.1, the velocity vo — 0.5 and = 0.01. Simulations 
(points) agree well with the theory (lines). 



decreases. In contrast to this behavior, Eq. ([^5]) shows 
that we recover for n = 2 a white Gaussian angle drive 
in the limit of small correlation times. 

For moderate correlation times a modified curled struc- 
ture emerges. For the stronger correlated systems one 
observes longer stays at a certain curvature 9 mixed with 
longer straightening segments where the random torque 
and the constant torque compensate each other. Thus we 
retrieve the pirouettes, but now with the preferred curl- 
ing orientation of the applied constant torques. Thereby 
the straight segments interrupted by the strongly curling 
pirouettes give rise to the maximal diffusive behavior and 
the resonance like structure shown with Fig. |6l 

Eventually, for very large correlation times we find a 
decrease of the effective diffusion coefficient, much like in 
the situation without torque [see Fig. ((S))]. DcH converges 
to zero for n — and to a non-vanishing constant for 
n — 1. The case n = 2 leads [i.e., due to Eq. ©] to 
a dominant torque in the limit of large Tc, so that the 
effective diffusion coefficient decreases toward 0. 

These asymptotic behaviors for small (i.e. ^Pq — ?► 0) 
and large correlation times (i.e. vI/q — >■ oo) in Fig. [5] can 
be discussed analytically in greater detail: According to 
Eqs. and ([25)) we arrive at the following limiting 

behaviors: 



if < f 



if *o > 1 



D. 



1 



off 



2D, 



« 1. 



D, 



cff 



(29) 
(30) 



Thus, we find a vanishing effective diffusion coefficient for 
correlation times approaching zero if n G {0, 1} and one 
converging to if n = 2. For large correlation times 

Tc, on the other hand, we find with 4' ^ 1 that -D^^ 



-15 
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FIG. 7: Spatial trajectories of an OUP driven angle with an 
additional angular force Q — 0.1, = 0.01, vq = 0.5 and 
with n = 0, for different correlation times Tc and the time 
length ti = 500; n = is used in order to obtain the largest 
Tc dependence of the displayed graphs (cf. Fig. |6|. 
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FIG. 8: (Color online) Effective diffusion coefficient DcS ver- 
sus correlation time Tc for different torques il at = 0.01, 
Vo = 0.5 and n — 0, within simulations (points) and in theory 
[eq. (O; lines]. 



converges toward zero if n = 2 or n = 0. In distinct 
contrast, Dcs approaches a constant value for n = I. 

In comparison to the asymptotic behavior without con- 
stant torque [cf. Eqs. ([25]) and ([?7)) ]. we notice additional 
multipliers with dependence on D, on the right side of 
Eqs. and (15111) . Both multipliers map onto the inter- 
val (0, 1). Thus, the effective diffusion coefficient always 
decreases with a constant torque V, ^ 0. This is depicted 
in Fig. |8] for n = 0. As both limits for large and small 
correlation times fall to zero, a maximal value for Dgg 
occurs, which not only decreases for growing but also 
shifts to larger Tc . 

We now inspect how noise and correlation time influ- 
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FIG. 9; (Color online) Noise intensity J)™^" which maximizes 
the effective diffusion constant for a given correlation time for 
different n at the velocity vq = 0.5 and the torque Q = 0.1. 
The inset shows the maximal Deft for these parameter pairs, 
which does not depend on n. 

ence each other for different n. We have seen that the 
effective diffusion constant exhibits a maximum for a fi- 
nite non-zero correlation time. Likewise, it has a max- 
imum for a finite non-zero noise intensity, as seen from 
the results with white noise (see Fig. [Ij. If we follow 
the ridge line of the maxima for fixed fl in the (tc, D^) 
parameter space, we obtain functional dependencies for 
different n-values, as depicted with Fig. [9l The inset 
shows the height at the ridge line (i.e., the value for the 
maximal effective diffusion coefficient as function of the 
correlation time). 

For small ft, we find that the optimal noise intensity 
jjmax^ which maximizes the effective diffusion coefficient 
at a given correlation time, behaves according to D^^'^^ cx 
t"~^. Further analysis shows that for increased 51 this 
relation gradually shifts to c>c t"~^~'^, with a small 

finite e. This shift can be justified by the denominator 
in Eq. (|2ip. in which the second contribution dominates 
for larger fl. Figure [9] depicts this dependence for an 
intermediate value of f2, since we observe a roughly linear, 
but still sub- linear, connection for n = 2, which converges 
to the value given for a white angle drive as Tc —> 0. For 
n = and n = 1, we recognize the expected reciprocal 
connection. 

If we follow the ridge line, at a given correlation time 
each has the same height for different n. We can under- 
stand this by considering that the parameters ^'o, t^Jc, and 
Dc, which contribute to the effective diffusion coefficient 
in Eq. (f20| . are all containing the same factor _D^t~" 
as sole n dependence [cf. Eqs. p3|) and ([TO])]. Hence, 
according to the discussion above, Dcs does not depend 
on n anymore while regarding D^^^. 

IV. CONCLUSIONS 

With this work the stochastic dynamics of active par- 
ticles with a constant speed that are additionally driven 



by an overall fluctuating torque has been studied. The 
cases of correlated and uncorrelated angle dynamics were 
considered and analytical results for the mean square dis- 
placement were derived. We discussed in this context the 
dependencies on characteristic parameters of our system, 
such as the correlation time Tc and the corresponding 
noise intensity D^, to find maxima in the mean square 
displacement. With respect to the correlation time and 
to the noise intensity a maximum of the effective diffusion 
coefficient D^fj was identified. We also gave a qualita- 
tive explanation of the maximal spreading of the agents 
by an inspection of the individual trajectories. The ran- 
dom composites of stretched and curved parts in the tra- 
jectories can give rise to an increased mean square dis- 
placement. 

It was also shown that the constant torque decreases 
the mean square displacement. Therefore, the persis- 
tence of motion and presence of a permanent torque 
are two appropriate instruments to optimize the two- 
dimensional motion of active agents. 

In view of the fact, that the displacement of an agent 
forms a strong link between theoretical and experimen- 
tal studies of active particles, we hope that further work 
within this field benefits from our theoretical reasoning 
presented here. 

Identifying the trajectories of our studied dynamics 
with the structure of a macromolecule, spin-off applica- 
tions in the field of polymer physics are conceivable. Es- 
pecially, the white Gaussian angle drive without a torque 
has its polymer analog in form of the famous worm-like 
chain model [3l|. Moreover, the similarities between the 
paths in Fig. [7] and Fig. 2] and correlated animal move- 
ments are striking for the case of animals, whose angle 
dynamics is determined by the past states. 

Finally, a useful extension of our model would be the 
adaptation to more realistic biological systems by ac- 
counting as well for an additional velocity dynamics v{t). 
Simple velocity models, wherein one can decouple the v- 
and 0-dynamics, are leading thereby directly to the dis- 
cussion in _3^, exhibiting a multi-crossover structure of 
the mean square displacement, which is due to the dif- 
ferent time scales within those systems. 
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